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ABSTRACT 

In the 't Hooft-Veltman dimensional regularization scheme it is 
necessary to introduce finite counterterms to satisfy chiral Ward 
identities. It is a non-trivial task to evaluate these counterterms 
even at two loops. 

We suggest the use of Wilsonian exact renormalization group tech- 
niques to reduce the computation of these counterterms to simple 
master integrals. 

We illustrate this method by a detailed study of a generic Yukawa 
model with massless fermions at two loops. 
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Introduction 



The dimensional regularization scheme devised by 't Hooft and Veltman jl], || 
and later systematized by Breitenlohner and Maison f| (BMHV) is the only 
dimensional regularization scheme which is known to be consistent in presence 
of 75 . It gives the correct result for the axial anomaly, at the price of breaking 
rf-dimcnsional Lorentz symmetry and chiral symmetry; while the former is easily 
recovered, it is a non-trivial task to satisfy the chiral Ward identities. 

More popular is the naive dimensional regularization scheme (NDR) M 
which, although inconsistent, is much easier to use and can be handled safely 
in most practical situations. For a review on the subject, see e.g. §| . 

Due to the relevance of higher loop computations in the standard model, 
where it is difficult to guarantee the consistency of the naive scheme, it is worth- 
while to investigate thoroughly consistent renormalization schemes. 

The difficulties encountered in computing systematically in the BMHV sche- 
me the non-invariant counterterms can be divided in three categories: 

i) it is a complicated (but generally solvable) algebraic problem to satisfy all 
the Ward identities (modulo anomaly problems), determining these countert- 
erms as combinations of Feynman integrals; 

ii) it is an analytically difficult problem to evaluate these counterterms, which 
in general involve Feynman integrals with several masses and/or momenta; ex- 
plicit formulae with several masses at more than one loop usually range from 
being quite complicated (see e.g. Q) to being as yet unknown; 

iii) it is burdensome to store these counterterms (which include evanescent 
ones) and to compute the renormalization group beta function using them j7j] . 

By comparison, we recall that in a vectorial theory (or in NDR, whenever 
possible) step (i) is trivial in the minimal subtraction (MS) scheme, since the 
Ward identities are automatically satisfied; step (ii) is much simpler, since the 
poles are much easier to compute than the finite parts of the Green functions; 
step (iii) is almost trivial, since in general multiplicative renormalization holds; 
simple formulae for the beta and gamma renormalization group functions are 
available in the MS scheme Qj. 

The short-cuts used to get the appropriate answer for the problem at hand 
include: 

I) in a few cases, like in the case of an axial current insertion in a vectorial 
gauge theory, it is sufficient to equate the axial vertex to the corresponding 
vector vertex multiplied by 75 to be sure that the chiral Ward identities are 
satisfied ||; 

II) after solving the algebraic problem, avoid carrying out step (ii); this is 
the philosophy of algebraic renormalization ( for a review see || ) , used generally 
only for demonstrative purposes, but which can also be used for making explicit 
computations, as shown in some two-loop examples carried out recently j JTo| . 

To our knowledge the only paper in which the counterterms are determined 
systematically at one loop in the BMHV scheme in a chiral gauge theory is jn|, 
in which the Bonneau-Zimmermann |L2[ identities are used. 

There are several computations of one-loop counterterms in the BMHV 
scheme for particular processes in the standard model JT3| , but no systema- 
tic one-loop treatment has been given in such a scheme. 

In this paper we describe a general method for simplifying step (ii) and partly 
(iii), i.e. we show that the counterterms can be reduced to zero- momentum 
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Fcynman integrals with the same auxiliary mass in all propagators, and that 
the beta function can be computed easily without making a direct use of the 
counterterms. 

The auxiliary mass technique in this form has been used in the past at two 
or more loop s fl4| , |l5| , |H| in conjunction with MS. Being the MS scheme mass- 
independent 1171 7 the auxiliary mass technique gives gauge-breaking terms that 
are polynomial in the auxiliary mass and then can be easily treated; however in 
chiral theories with BMHV the MS scheme cannot be used, and it is not easy 
to disentangle the gauge-breaking terms. 

We will show how to do this by Wilsonian methods fllgfl , 
Wilsonian methods have been used by Polchinski Jl9| to simplify the proof of 
renormalizability in <f> A \ this proof has been further simplified and generalized to 
other cases ||(], ; in particular gauge theories have been renormalized 

using the effective Ward identities introduced in [^(J. In mass- independent 
renormalization has been studied with these Wilsonian methods; it has been 
also shown in this paper that the Wilsonian effective action satisfies an effective 
renormalization group equation, which is the analogous of the effective Ward 
identities. 

In this paper we propose the exploitation in the BMHV scheme of the effec- 
tive Ward identities and effective renormalization group equation to compute 
the finite counterterms and the beta function in terms of Wilsonian Green func- 
tions at zero momenta and masses, which are easy to evaluate. 

As a simple testing ground for our proposal we renormalize systematically 
at two loops the most general Yukawa theory with massless Dirac fermions; 
so far only the simplest Yukawa model with pseudoscalar coupling (without 
chiral symmetry) has been renormalized at two loops in the minimal BMHV 
subtraction scheme . 

We impose Wilsonian mass-independent renormalization conditions compat- 
ible with the rigid chiral Ward identities. Our renormalization scheme is chosen 
to coincide with the minimal subtraction scheme in the non-chiral Yukawa case 
(this choice is done to simplify step (i), but is not essential in our method); in 
the general case it gives the same two-loop f3 and 7 functions as in the minimal 
naive scheme |p4) . 

Finally we discuss the coupling with external gauge fields; we impose Wilso- 
nian renormalization conditions compatible with the effective Ward identities 
at one loop. At two loop we check only the gauge two-point function Ward 
identity and the Adler-Bardeen non-renormalization theorem |25| . 

Our method is based on the use of a Wilsonian flow belonging to the class 
characterized by the cut-off functions K^ l \x) = ( X + A z ) n with n = 2,3,...; 
these cut-off functions separate the propagators into hard and soft parts. For 
n > 2 the counterterms chosen at a particular A renormalize the flow for any A. 

We renormalize the theory at zero momenta and masses along the flow n = 2, 
using the n — 3 case as a check on computations. 

The counterterms are zero-momentum integrals with the same A in all the 
propagators, which can be reduced to a single 'master integral' at two loops 
using recursion relations (for a review see |l5| ). We have used Mathematica 
to evaluate these integrals. 

In the first section we show how the exact renormalization group method 
can be used in connection with dimensional renormalization to renormalize the 
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Yukawa models. 

In the second section, after reviewing the recursion formula for the massive 
zero-momentum two-loop integrals, we give the one and two loops results for 
the Yukawa models. 

In the third section we discuss the effective Ward identities for the Yukawa 
model in presence of external gauge fields. 

We end with a concluding section. 

In the appendix we give the complete two-loop bare Yukawa action. 

1 Dimensional renormalization and the Wilso- 
nian effective action 

1.1 't Hooft— Veltman regularization scheme 

We recall how gamma matrices are treated in the 't Hooft and Veltman jl] 
dimensional regularization scheme as elaborated by Breitenlohner and Maison 
1 (BMHV). 

We work in Euclidean space; <5 M „ is the Kronecker delta in d = 4 — e dimen- 
sions; 

5p V p v =Pp ; S l±l± = d (1) 
The gamma matrices 7 M satisfy the relation 

{y^yu} = -26^1 (2) 

where 

trI = A ; 7 7 „= 7jU / = 7 „ (3) 

In the BMHV scheme O(d) invariance is broken and one introduces 0(4) x 
0(d — 4) invariant tensors : the (d — 4)-dimensional Kronecker delta 5^ and 
the 4-dimensional antisymmetric tensor t vp(J T , satisfying 

^pv^vp , ^pp 6 , 6 ^u^u par (4) 

Moreover one defines the 'evanescent' tensors 

Pp. = ; > ee (5) 

The matrix 75 is defined by 

75 = jfpvpalplvlpla (6) 

satisfying 

{75,7mI = {75,%} = 2%75 ; 75 = I (7) 

In Euclidean space the reflection symmetry takes the place of hermiticity. 
Reflection symmetry is an antilinear involution, under which 

®1p = -1p ; 9^(x) = $(a/)>n ; ®${x)=^{x') (8) 
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where x 1 = — 



x 1 , x = for fj, ^ 1. 



For a multiplet of fcrmionic fields, a general local marginal fermionic bilinear, 
involving a scalar (A), pseudoscalar (B), vector and pseudovector (A^), all 
of them real fields, can be written as 



where Hi are matrices over flavour indices. Reflection invariance requires that 
Hi are hermitian. 

Green's functions in d dimensions are obtained performing the Dirac algebra 
in the Feynman graphs with the above rules. Actually the analytic continuation 
to continuous dimensions is defined on the scalar coefficients of the Green's 
functions expanded on a basis for tensorial structures. For chiral theories, being 
0(d) broken, such a basis includes also 'hatted' tensors. 

The poles of the Green's functions for e — > are removed by local countert- 
erms. Loop by loop the singular part of the counterterms must subtract exactly 
the pole part of the Green's functions, including the 'hatted' components; the 
finite parts of the counterterms instead are not uniquely determined by the 
renormalization conditions which constrain only their 4-dimensional part. In 
fact the d — > 4 prescription requires, besides to take the limit e — ► 0, to set the 
hatted tensors to zero. 

In order to fix completely the fcrmionic counterterms, in the decomposition 
(ph we will choose that H3, H4, Hg, H\q have vanishing finite term in their 
Laurent expansion. The same convention is assumed for the coefficients bij, dfj 
and df!? in the bilinear scalar counterterms, which have the form 




4>{Hi-f^ + Ha[7/i) 75] + #37m + «#47m75)cV0 + 



ip{iH h A + H Qlb B + iHrypVp + i-ffsfo-pi 7s] ^ + 



(9) 



iH^vVfi + #107^75^)^ 




(10) 



1.2 Yukawa models 



Consider the most general Yukawa model with massless fermions. 
In dimensional regularization the bare action is 




(11) 



and it is chosen to be reflection symmetric. 
At tree level 



c ij (p) = 6 *jp 2 + m % ; ' 

(0) _ ■ . (0) _ 7 
c i — l Vi i c ijkl ~ n ijkl 



c {0 '(p) = z7 M p M 



(12) 



We define 



yi = Sil + iPi7 5 ; Yi = Si + iP l 



(13) 
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The constants c, Cj are matrix- valued, c = c/j , etc., where i,j are the internal 
indices of the fermions ipi- The matrices Si, Pi are hermitian. 

We will consider a group G which is not necessarily semi-simple, with struc- 
ture constants f abc . The fields transform under linear (chiral) representations 
of the group which are in general reducible: 



where 



and 



if) — >• gtp ; V -> V>5 1 ; -* ^i<^j ( 14 ) 
gr = ejrp(ie t°) ; g = exp(ie a P) ; h = exp(ie a 6 a ) (15) 



t a = t R P R + t a L P L =t a s +t a pl5 ; t a =t a L P R + t R P L =t a s -t a pl5 

1 1 (16) 

ft =2(1 + 75) ; ft = 3(1-75) 

t R and i£ belong in general to different representations of G. The scalar fields 
are in a real representation; then 9 a are antisymmetric imaginary matrices. 

The Yukawa coupling ipyiip(f>i is invariant under these transformations pro- 
vided the following relations hold 

y j e° i + y i t a -t a y i = ; y)d% + y\t a - t a y\ = (17) 

or equivalently 

y 3 0%+Y t t R -t L Y^Q ; Yje^ + Yjtl-tpr? =0 (18) 

The tree-level action is invariant under these chiral transformations, apart 
from an evanescent fermionic kinetic term. Higher corrections in the bare ac- 
tion will require also non-invariant counterterms; in the next subsection we will 
discuss how Wilsonian methods can be applied to determine the counterterms 
in order to preserve the Ward identities. 

1.3 Wilsonian effective action 

The bare action (|ll|) has the form 

= + S 1 ^) (19) 

where we use a compact notation in which $ is a collection of fields. S 1 contains 
the tree-level interaction terms and the counterterms. 
Consider an analytic cut-off function 

for n > 2 ; M 2 = m 2 > in the scalar sector; M 2 = in the fermionic sector. 

Let us split the propagator in two parts characterized by a scale A > 0; 
defining 

D H = D A = D(l - K A ) ■ D S = DK A (21) 
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Let us make an 'incomplete integration' over the hard modes 



Z A [J] = exp~W&.[J\ = J 2?$ e -^ 5A W- J *l (22) 

with bare action 

S A ($) = ^DX 1 $ + S I ($) (23) 

The Green functions obtained from Z A [J] are infrared finite for A > even at 
exceptional momenta. 

The flow of this functional from A to zero can be represented as 

Z\J\ = Z a [J] = e^^ D ^ KK ^ Z A [J] (24) 

The Wilsonian theory at scale A has a bare action differing from the one of 
the usual theory at A = by the term 

5a(«) - = ifcDX 1 ^* (25) 

which has ultraviolet dimension less or equal to zero, due to the above made 
choice of the cut-off function K A ; the renormalization of the usual theory implies 
the renormalization of the Wilsonian theory and viceversa. 

To impose the renormalization conditions it is useful to separate the Wilso- 
nian 1-PI functional generator, obtained by making a Legendre transformation 
on VFa[J], into the quadratic tree-level and interacting parts: 

r A [$] = i$^ 1 $ + ri[<i>] (26) 

and to make a Taylor expansion of r^[<f>] in fields, <i-dimensional momenta and 
masses; the four-dimensional terms in this expansion form a local functional, 
which will be called Sw- 

Let us now discuss briefly why the renormalization program can be per- 
formed on the Z A functional with some advantages: 

i) since the hard propagator Dh (p) is regular at p = there are no infrared 
problems in the evaluation of the counterterms, so that the counterterms can 
be chosen to be Feynman integrals at zero momenta and masses using a mass- 
independent Wilsonian renormalization scheme p3| ; 

ii) for a Wilsonian flow with n > 2 it is simple to prove that the renormaliza- 
tion at a fixed A implies the finiteness of the Wilsonian effective action at every 
value A of the Wilsonian flow, in particular at A = 0; moreover because the 
propagators D\ at A = do not depend on n, a bare action that renormalizes 
a given flow will also make finite the flows for every n > 2. 

iii) for a Wilsonian flow with n > 2 the Ward identities and the renormaliza- 
tion group equation on the functional Z A= $ are equivalent to the corresponding 
effective Ward identities ^2) and effective renormalization group equation 
[ p3[ on Z A , so that choosing renormalization conditions compatible with the 
effective Ward identities the validity of the Ward identities on Z A=0 follows; 
furthermore the renormalization group beta function can be expressed in terms 
of Feynman integrals at zero momenta and masses. 
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These considerations hold for any renormalizable theory in four dimensions; 
we consider here only the Yukawa model, which is the simplest chiral theory. 
Actually its rigid Ward identities are too simple to illustrate point (iii), since Z\ 
satisfies the same rigid Ward identities; however introducing external currents 
one has to study the local Ward identities, whose effective form will be used in 
a later section. 

In the Yukawa model the bare constants in the action are chosen of the form 

(27) 

r>0 

We use a modified subtraction scheme [^7j , in which Nd is introduced in 
order to avoid 7^ , ln(4ir) and tt 2 factors. We choose a mass-independent 
renormalization scheme, i.e. all the bare terms, apart from Cy, are independent 
from m 2 ; the latter term depends polynomially on it. 

The bare constants must be fixed by suitable renormalization conditions on 
the Wilsonian effective action. 

In the Yukawa model the marginal part ^ of the Wilsonian effective 
action is defined according to: 



t^A rtA 1 tiA 

1 — >~>W "T 1 irr 

/- 1 1 

^ad^ip - -aijcfiid 2 ^ + -(^to^X/^ + (28) 

1 

IT 

Ti rr contains all the terms of dimension greater than four in the expansion of 
Ta in the fields, momenta and masses. 

The renormalization conditions fix the limit for e — > of the functions — 
{a, ciij ;, a^J, cii, aijki} at A = [i and perturbatively determine uniquely the bare 
action, as discussed in the first section. 

From the Feynman graphs rules the terms gA at order I in the loop expansion 

are 



.9,1 



^A' 

r=0 



U) =E(T) r ^ (29) 

in which the dependence on A and /i is made explicit; g^h IS the '-loop graph 
contribution, while g^' 1 ^ , r = 0, I — 1 is the contribution due to the r-loop 

graphs with counterterms of overall loop order I — r; g^' 1 r ^ are independent 
from fi and A . 

Being the theory renormalized, one gets, for e — > 



^ = - A ^^ r ti" r) (30) 



r=l 



2 In a Wilsonian mass-independent scheme, as discussed in |23| , the renormalization of the 
mass parameter is treated in a way very similar to the kinetic term and therefore is included 
in the marginal part of the action. 
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Taking s derivatives with respect to A in (|29|) one has still a finite expression 
for e — > 0, so that one obtains the consistency conditions 



Xy^^o (si) 

r=l 

for n = 2, I and s = 1, n — 1. 

For the two-loop case the consistency condition ([H]) for n = 2 and s = 1 is 
known to hold for each Feynman graph and its counterdiagrams p8| , providing 
a useful check. 

To renormalize the theory we assign finite values to the constants qa at 
A = /i 

9A = r A + 0(e) (32) 

Each coefficient in the Laurent expansion in e must be determined to estab- 
lish completely the renormalization scheme. In a subtraction scheme in which 
the bare constants are fixed to have all positive powers of e equal to zero (see 



eq. (27)) and the renormalized quantities are fixed to be non singular and with 
determined e° coefficient, the renormalization scheme is completely fixed. The 
0(e) term in eq. (|3^) refers to the fact that one cannot fix the positive powers 
of e both in the renormalization condition and in the bare action. 

In order to respect the rigid Ward identities one must choose for ta group- 
covariant quantities under the chiral transformations. They can be constructed 
out of the tree coupling constants y and hijki, taking into account the invariance 
properties 

g' 1 Vig = hijyj ; g~ 1 yjg = h^y] (33) 

Furthermore one can examine which choice for ta gives the simplest expression 
for the countcrterms. 



For the non-chiral theory, in which the couplings Pi of eq.(13) are vanishing 
and g — g — g s , the simplest counterterms are those determined by the MS 
scheme, corresponding to renormalization conditions ta{S) that can be explic- 
itly computed. 

For the chiral theory we will define the counterterms as suitable functions 
ca(S, P) which, for P = 0, coincide with the corresponding functions of the non 
chiral case. To this aim loop by loop we will choose ta applying a covarianti- 
zation formula to the corresponding ta(S). A comparison of eq.(^3|) with the 
analogous equation in the non chiral case with the same group G: 

g~ 1 S i g s = h ij S j (34) 

suggests the recipe for the Green's functions of interest : 

Si 1 Si 2 . . . Si 2n + 1 I > yiiUi„ ■ ■ ■ Vtm+i 



jf,S tl S i2 . . . S l2n -> 7^*4 y i2 . . . y i2n 
1 
4 



Tr [S il S i2 ...S i2n ] -» -Trtr 



Vi\Vi 2 ■ ■ ■ v\ 2 „ 



(35) 



1. 



-TV [Y l ^..Y l2n _^l n +Y}Y l2 ..Yl n _Y l2n 



where Tr denotes the trace over the internal fermionic indices. 

As an example let us consider the renormalization condition to be imposed 
on the fermionic self-energy. In the non chiral theory, using the minimal sub- 
traction one computes: 

S(p)| roars =i$r(S) (36) 

Using the covariantization procedure ip r(S) must be replaced by ij> r(y, y'). In 
order to be compared with the standard form of eq.(|J) S can be written as 

Z(p)\rnar a =^U [v(Y^) + r(Y\Y)} + 

21 (37) 
-0,7b] [r{Y^)-r{Y\Y)}} 

which differs from vp r(y, y^) by evanescent terms. Eq.(|37]) gives the renormal- 
ization condition for the chiral theory. Observe that, if the theory is reflection 
symmetric, r(y,y' ! ) is hermitian. 

The counterterms ca are completely determined by the pole part and by the 
constant part of the corresponding vertices. They can be decomposed into two 
parts 

ca = c N A DR + Ac a (38) 

where c^ DR has the same group structure as r A and, due to the choice of 
renormalization conditions, has only pole part in e (see comment after eq.(|32|)). 

The remaining part Ac A vanishes in the non-chiral case (Pi = 0). 

The non-marginal relevant terms satisfy 

r^ _0 | p=m= o = ; rfj \ p=m=0 = (39) 

In dimensional regularization, being the massless tadpoles equal to zero, 
eq.(|39|) corresponds to have vanishing bare relevant counterterms. 

2 Explicit computations 

2.1 Master integrals 

Using the cut-off function ( pp| ) the hard propagator has the form 

D A ( P ,M)=D( P ,M) P -^Y: L 2 + M2+A2 ) (4°) 

r— 1 ^ ^ ' 

The counterterms are computed in terms of Wilsonian Green functions at zero 
momenta and masses M — 0; the corresponding Wilsonian Feynman graphs 
with / internal lines have I sums 2™ =1 •■■ E" ; =i- 

Since D (p, M)(p 2 +M 2 ) is polynomial in p the counterterms are expressed in 
terms of massive zero momentum tensor integrals, in which all the propagators 
have 'mass' A. Taking traces one can reduce these tensor integrals in terms of 
scalar integrals. 
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la = I -TT^T^r-^ (41) 



At one loop the scalar integrals are 

d d q x 1 

At two loops the scalar integrals have the form 

_ f Ai f d% 1 

aAc J N d w^ J N d <K^ (qf + lY{ql + l)>>(( qi + g 2 ) 2 + l) c 1 ' 

and so on at higher loops. 

At one loop one has the well-known explicit expression 

N d T{a) 

At higher loops an exact expression for these integrals is not known; they 
can be reduced, using recursion relations, to a small number of master integrals 
which can be expanded in e; to renormalize at Z-loop one must know the 1^ 
up to the O(l) term, the up to the O(e) term and so on (in minimal 

subtraction it is sufficient to know the 1^ up to the 0(l/e) term, the /C -1 ) up 
to the O(l) term and so on; however in presence of chiral symmetries MS is not 
sufficient). 

At two loops there is only one master integral (for a review see e.g. flStO ; it 
is known how to compute the finite parts of the three-loop master integral ||16|] . 

In this paper we will restrict to two-loop computations; the recursion relation 
is obtained from the identity 

d d qi f d d q 2 d gg _ Q 



N d lT d /2 J N d ^ dqt {q\ + l) a {q 2 2 + l) b ((<?i + q 2 ) 2 + 1) 
which implies 

(d— 26 — c)I a +l,b,c — c(I a ,b-l,c+l — Ia-l,b.c+l) + 
26/o,6+l,c + cl a ,b,c+l = 

From this relation and the fact that I a ,b,c is totally symmetric in its indices it 
follows the recursion relation 

3ffl/ a +l,6,c = c (Ia-l,b.c+l — Ia,b-l,c+l) + 

b {Ia-l,b+\,c — Ia,b+\,c-l) + (3a ~ rf)^a,fc,c 

Using this recursion relation one can express all I a .b.c , for a, 6, c > in terms 
of the master integral 71,1,1 and in terms of I a >,b'.c', with one of the indices 
vanishing; in the latter case the two-loop integral reduces to the product of two 
one-loop integrals: 



I a ,b,o = Ub (47) 



One has 



fi 9 3 

h,i.i = -4-- + ^- 5 )+°( £ ) 



4 r /3 / x\ (48) 

v = 3 lim I 2 2 2 = -2 = / dx In ( 2sin- ~ 0.34391 

e^O ' ' 73 Jo V 2 J 

These recursion relations can be solved very fast by computer in the cases of 
interest. Recently all the Laurent series of /i,i,i has been computed in J29|. 
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Table 1: One- loop coefficients 



2.2 One-loop results 



At one loop the renormalization conditions (32) on the Wilsonian vertices are 



n h ■ ■ 
(l) "y 



J 1 ) _ r <t> Y- ■ ■ r^ rs — *- m 

'ij ~ '(I) 1 *! ' 'ij 

(1) if> t (!) • v3 t 



(49) 



r ijfei - r (l)i r (yfcO + r (l)2 n (ij n . 



kl) 



where the symmetrizations (...) are with weight one: e.g. h^j k ^ = hijki; we 
defined 



Y 



1 



Hl2~ '2n-l*2« 



= -Tr V, V/...K 



*1 J ? 2 ' ± l 2n-l 1 i 2n ~ il l 2" «2n-l * 2 



(50) 



The same parametrization holds for the bare constants cy[ defined in 

The quantities ta depend in general on the choice of the Wilsonian flow; in 
Table |l] we give the values of ta and Ca DR for the flow n = 2, figure [l] gives 
the corresponding graphs. To compute the coefficients in Table yj one has to 
make an expansion in masses and momenta of graphs in fig.jj], as mentioned 
after eq.@ 

The remaining part of the one-loop bare constants is 



A 4) } = i ( -p 2 + -f ) Tr [PiPj] 



Ac« = -j> (2P i P i -ij 5 {S i ,P i }) 



Ac« = 



A 4]li = -96 Tr 



S^PkPi) + ^P(iPjPkPi) 



(51) 
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(e) (f) 

Figure 1: One- loop graphs 

In 0] MS is applied in BMHV in the simplest Yukawa model with a pseu- 
doscalar, where there is no chiral symmetry to be maintained. Here we intro- 
duced finite counterterms, which are exactly those needed to obtain the same 
renormalized Green functions as in the MS NDR scheme. In it was shown 
that the beta function at two loop in the MS BMHV scheme differs from the 
one in MS NDR scheme by a renormalization group transformation involving 
finite one- loop counterterms, which are in agreement with eq. (|5l|) . 

2.3 Two-loop results 

At two loops the renormalization conditions on the Wilsonian effective action 
read 

d> 2 i m 2 ,kl 2 
r v = r ijP + r n ' m kl 



13 

m 2 ,kl 



(52) 



— r m h , h -, 4- r m h ■ 
' 1 ' L ikmn' L jlmn ~ ' 2 i. 

f^ 1 Yikji + f™ Yijki + r™ Yi m hijkm 
The two-point fermionic term is 

r = v)\ r tyiv\vi + rtviVjyi + rfYij y t + rfh^kkVi] (53) 
The two-loop fermion-fermion-scalar term is 

Ti =iyk[rfY jk ylvj + rfyjyiyjyk + rf{yjyjy] + y]yjy})y k + 

rf(y\y 3 y{ + v 3 WJ)l/j + rfy]y^yj + rfh ]k uy\y d + rfh tjkl y]yi] 

The two-loop quartic scalar term is 



(54) 



Tijkl — J-ninjkl ~r 7*2 *nijnkl ~r 7*3 1-nijkln \ 

hmnij ^rnkln H~~ ^5 ^mknl ^rnphpnkl H~ (^^) 

hmri'pqh"pqkl H~ 7*g hmpklhnpqq H~~ 7*9 ^npqk hmpql) 

where symmetrization over the indices A*, I is understood. 
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Table 2: Two-loop coefficients for the two-point bosonic functions 
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Figure 2: Two-loop graphs for the two-point bosonic function 
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Table 3: Two- loop coefficients for the fermionic two-point functions 
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Figure 3: Two-loop graphs for the two-point fermionic function 
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Table 4: Two-loop coefficients for the cubic vertex 
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Figure 4: Two-loop graphs for the cubic vertex 
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Table 5: Two- loop coefficients for the quartic vertex 
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Figure 5: Two-loop graphs for the quartic vertex 
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The coefficients in eqs.(|52j) , (153) , (M), (|55j) and the corresponding ones for 
the two-loop c^ DR in the case of the Wilsonian flow n = 2 are given respectively 
in the Tables g, ||, [|, |[ Figures ^, [| Q, show the graphs contributing to these 
coefficients; for each diagram shown there are in general corresponding graphs 
with one-loop counterterms that are not represented. 

The naive part c^ DR of the two- loop bare constants agree with (24j . In order 
to make the comparison with VM, where a two-component spinor formalism 



is used, one must substitute the traces kTr 
Y 

1 Ill2-..l2n-ll2n ' 



Y Y Y^ 



with 



The remaining part of the bare counterterms is given in the Appendix. 
With the bare action determined using the flow n = 2, according to the 



lines of sect. (1.3) one can define a new Wilsonian action T^ 1 . For n = 3 



(or bigger) the new effective action is finite and satisfies the Ward identities. 
The graphs of T^ 1-3 ' now have n — 3 propagators D\ but contain the n = 
2 counterterms. This suggests the following consistency check: consider the 
contributions of n = 3 bare (marginal) graph, of its relative countergraph and 
two loop counterterm; we have verified that although the two first quantities 
are divergent and separately not invariant, the total sum is finite and invariant. 
Furthermore this total sum, being computed on a different flow is different from 
the n = 2 case. 

2.4 Renormalization group equation 

Let us compute the renormalization group beta and gamma functions in the 
BMHV scheme at two loops. In the case of the simplest Yukawa model this 
has been done in Q] using the minimal subtraction formulas || and taking into 
account the one-loop evanescent tensors in the bare action; after obtaining a 
renormalization group equation involving also insertions of evanescent tensors, 
these are solved in terms of relevant couplings obtaining the usual renormaliza- 
tion group equation. 

To avoid making a similar subtle analysis of the bare couplings in our case, we 
will obtain the beta and gamma functions working directly with the Wilsonian 
effective action, which satisfies the effective renormalization group equation p3|] . 
We will restrict to the massless case; in [^3| it is considered also the massive 
case. 

In the compact notation of subsection 1.3 the Gell-Mann and Low renormal- 
ization group equation reads 

<+^| + ^W]=^] (56) 

where Z is the renormalized functional in which the limit e — > has not yet been 
taken; £ is an evanescent functional, which in general is not vanishing for finite 
value of e since the renormalization of the theory is not strictly multiplicative, 
and 

"■5-*»s4r +lvw,i ^ + *^ ) (57) 

From eq.(E3) it follows that Z\[J] satisfies the 'effective renormalization 
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group equation' 

Define the functional Z\[J,x] as in (p2|), but with 

5 A ($,x)=5 A (<i>)-xA 7 [<i>] ; A^-^D^K^ (59) 
The effective renormalization group equation can be written as 



(v—.+P- — - - — )| x =o)r A [$,x] = £ A ffl (61) 



which in terms of the 1-PI functional generator reads 
or equivalently 

{ ^ +P '§-g- ^JQ^W = + Em (62) 

where T^^] represents the insertion of the non local operator A 7 [<1>] on the 
functional r A [<!>]; its Feynman graphs contain the new vertices 

■ 2K A (p)D~' L (p) 7ij 

■ ^ 2K K {p)S H 1 {p) 1 ^ 

between two bosonic or fermionic lines; due to its nonlocality A 7 [$] does not 
require renormalization. 

Let us finally rewrite in a suitable way the gradient terms in (p2|); using ( |30| ) 
one gets 



These gradients have the same group structure as the renormalization constants 
given in eqs.(E9|), (p2[); the coefficients of these group structures are given in 
the tables [l], [§ U § hi the case of the n = 2 flow. Observe that, while 
<?a'.?i an d <?a'_i are n °t separately chiral group-covariant, the gradient term 

2<?x — l "I" 5a —l mus t be covariant. This provides a check on the computations, 
besides the double pole rule mentioned after (31), that we made in a systematic 
way. 

At one loop the Wilsonian gradients can be expressed in terms of the pole of 
the bare couplings, so that the formulae for the beta and gamma functions are 
equivalent to the standard formulae 0; using Pa — (ifp $a an< ^ defining 

Pi — Py z Pr + 0y. Pl we get 

J 1 ) - 2Y ■ - -7/-i/ t • y^ 1} - -vUj- 

= IVjViVi + \(ViVjVj + yjVjVi) + 2y y% (65) 
4-1 = -48y (yw) + 3^?-^ + 8Y m(t h m m 
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Figure 6: 4-point T Green function contributing to Ynijkin 

At two loops the Wilsonian gradients are not given by the simple pole of 
the bare counterterms (|64|); the beta and gamma functions are given in terms 
of these Wilsonian gradients, the T insertions, given in the tables § |, § [| and 
the one-loop renormalization conditions (f49|). 

We get 



(2) 
lij 


— -. ~hikmnhjkm,n 
12 


^^ikjk 


^^ijkk 




1 t t 

= -^jVivlyj - 






7 ^(2 


1 t t 

= —gVMViVj - 






(IP 


= 2y k y]yi{-y]yk 


+ vlvj)- 


Zhijkiyjylyi - 4Yj k yjy}y k 



(66) 



ykiViVjVj + y)yjy\)yk + yn i ' {2) + 7* (2) ^ + i\f 



(2) 

$ijkl — 96(21 / ^ l 2 n j/ c / -|- Yfiijnkl H~ ^nijkhi) ~\~ 



(2) 

where in the last equation the indices fc, i must be totally symmetrized. 

These formulae are in agreement with p4j . No renormalization group trans- 
formation is needed, because the renormalization conditions are the same. 

As an example, we give the separate contributions to the term Y n ijkin of 
Pijli ■ I n (H) are indicated the Feynman rules for the insertions . 



-/i-fjf 1 (see table |) contributes igjp 



• ^9ij k i gets the contribution -^f^ due to the term ^ Y(ijH) of 9ijli 

(see table 0) and the term \{yiy]yj + Vjyjvi) of /3f\ 

Finally 7^ contributes to Ynijkin through the graph in fig. || and gives the 
coefficient "y^- 

These three contributions add up to the expected value 96 Y n ijkin for p\^ kl 
in equation 



3 External currents and local Ward identities 
3.1 External currents 

In order to define the currents associated to the symmetry transformations (jl4j) 
we will add in the action ( jll| ) an external gauge field A". 
At tree level the action is now 

f 1 - 1 

5 (0) = / ipi^D^if) + — (D a <j)) 2 + iipyiif>(f>i + ^jhijfdfctjijfafa (67) 
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where 



D^ = (d^ + iA«t a )i; ; = {d^ + iAffity (68) 

We can choose for the external gauge field the condition A a ^ — 0, so that 
Jnt a is equivalent to i a "f^ and the fermionic gauge coupling «7 M i a is equal to the 
reflection-symmetric expression |(7 p i a + i a j^). 

We set the // e factors to one; they can be easily reintroduced by dimensional 
analysis. 

The local infinitesimal version of transformations dl"i] ) is now 

8 A* = -d„e a - f abc A c ^ b = ~D ab e b 

Sip = i( a t a i> ; Sip = -ie a 4)i a (69) 

Stk = 'V'W'; ; o, 

where d^e a = 0. The tree- level action is invariant under these transformations 
apart from an evanescent fermionic kinetic term. 

To renormalize arbitrary products of currents we will consider in the bare 
action local monomials in A a ^ up to fourth order. 

The bare action (O) gets additional terms: 



(70) 



/1 1 1 
-r ab A a A b -l r abc A a A b A c -l r abcd A a A b A c A d 4. 

AtfW + -c^&fc + \c;1 3 a;a%^ 

The marginal part Sw of the Wilsonian effective action will have a similar 
dependence on A", 

The renormalization conditions on Sw analogous to (|3^) , which include the 
new vertices, have to be chosen in such a way that the local chiral Ward identities 
hold; to this aim the renormalization conditions are fixed compatibly with the 
effective Ward identities. These identities have been introduced in ^(J and 
further studied in fE2|. 



An analysis similar to the one in sect. 2.4 shows that the extra term of the 
effective Ward identity can be represented as the insertion of non-local vertices. 
Under transformation (p9j) one gets: 



ST A =T A +(D A 7~ = 



^$D^K A <Z> 



• r A (71) 



where T A is the generator of the 1-PI Green functions with one insertion of the 
operator S (i<f>Z? j ^ 1 i^A^ > ) • Its Feynman rules are 

- ie a (p)[S^(p + q)K A {p + q)t a - t a (q)K A (q)} (72) 

for insertion on the fermionic lines and 

- ie a (p) [Dg 1 (p + q)K A (p + q)6 a - e a D H x (q)K A (q)} (73) 

for insertion on the scalar lines. 

If one can choose the renormalization conditions on T A such that O a is 
evanescent, the effective Ward identities are not anomalous and the usual Ward 
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identities are recovered at A = 0; otherwise it is anyway possible to separate 
the functional O a in a part corresponding to the local anomaly operator plus a 
part representing an evanescent operator: 



(74) 



cast in a form which we 



Equation (j7l]) is the effective Ward identity 20, 
find convenient for explicit computations. 

The anomaly operator is defined by its renormalization conditions at the 
rcnormalization scale A = fi therefore we have to consider the marginal projec- 
tion of eq.(|7l"l) in the limit e — *■ 0: 



SS- 



w 



T + A 



(75) 



where T and A are the marginal parts of T A and A A . 

The additional rcnormalization conditions to which we refer in the following 
sections, have been chosen through a covariantization procedure of the renormal- 
ization conditions of the corresponding vector theory with minimal subtraction. 
For example to the tensor Tr[t^Sit b Sj] we associate in the chiral theory the ten- 
sor jTrtr[t a yli b yj], which has the same transformation property. According to 
this rule by construction Sw does not contain monomials in the gauge fields 
with the Levi-Civita tensor. This leads to the definition of the anomaly in the 
left-right symmetric form. 

The Bardeen anomaly can be obtained through different renormalization 
condition, allowing in Sw terms proportional to the Levi-Civita tensor. 



3.2 One loop 

At one loop using the tree level action ([ 
one computes: 



and the rules (|72|,|73|) for the insertions 



T 



J V e ^>?Dfd v Al + bfr b Al{d^A b v - d v Al) 

(febx b xyjyac + ^abc^a^c + 



A 



(76) 



-4 = / -ie^d^Tr 



t e R (d v A Rp A Ra + -A Rv A Rp A Ra ) - (R -> L)J (77) 



where A Rfl = it a R A«, A Lp = it a L A<^. 
With the notations: 



1, 



S 2 (Fr = -Tr[t R t R + t a L t b L ] 



S 2 (S) 



ab na nb 

- <>,.>>,. 



\abcd na nb nc nd , na nb nd nc , na nc nb nd 



S 4 (S)° 



(78) 



20 



the values of the coefficients in eq.([76|) are: 
b? = -ls 2 (F) ab + ±S 2 (Sr b 

bf = -ls 2 (F) ab + ±s 2 (sy b 

34 7 

b abcd = AS ^ F ybcd _ }}_S 4 (S) abcd 

„ 7 10 (79) 

b % = J Tr K y /*L + Y}YjV& - (* ^ i) 

1 Oj , 

6 « = 210 ( Tr K*&'*k> y < ty i + + (*" 3) 

i| (rr^r/^y, + t^r/] + (» - i) 

To determine the Wilsonian renormalization conditions one must find ; 
choose them of the form 

S [ w ] = JlafA^Al + \a a 2 h Ald^ v A b v + a a 3 x f^A^Al + 

h f ac Xa xv f ybd + l a abcd )A a A ^ A c A d + iA ^ a ^ + (80) 



6 

1 

4 



The effective Ward identities give the following relations between the S^P 
coefficients and those for T : 

(oi + a 2 + bi) ab = (ai + a 3 - b A ) ab = 

(03 - a 4 - b 6 ) afc - (a + 6) afccd = 

t e a-a e + b e = [i e ,a]=0 (81) 

(a" 4- 6 a - ie°o)« = ia? k 6 b kj + ia%0 b ki - a$ ~ b b f = 

Using (|8l|), a set of independent parameters of Sy(P is a, ay and a^ 6 . The 
former two parameters have been fixed in the previous sections, to be compatible 
with minimal subtraction in the non-chiral case; making the analogous choice 
for the new parameter we get 

af = \s 2 {F) ab + ^-S 2 (S) ab + z ab (82) 

with z ab = for the minimal subtraction choice. We have checked all the 
relations ( |8l| ) by explicit computation of the marginal Wilsonian Green functions 
and the corresponding T insertions. 

Having determined the renormalization conditions, we can now compute the 
one-loop bare action: 



S(i) = J \F^[~S 2 (Fr b ~ ^ 2 (5) Qb + z ab ]F b u + 



G cj>NDR . G ipNDR . a a 
b tl) + b (l) + AA (1) 



(83) 
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where F° u = d^A a v - d u A% ~ f abc A^A c v , Sf^ DR and Sff DR are the covari- 
antization of the kinetic terms for fermions and scalars in absence of the gauge 
field; AfffD is the non-naive part of "S'(i). 

A in ([77]) is the Adler-Bardeen anomaly (2^] ; for a previous derivation using 
Wilsonian methods, see p^j. 

In order to compute the Bardeen anomaly in the next subsection we write 
A5^\ the gauge-dependent part of AS(i), with the variables and Vp so 
defined: 

A„ = V ll + 'y 6 A lt ; V„ = it a s A^ ; A^it^A" (84) 
As a consequence: 

= - d v Ap + [A^, A v ] = Vp V + j 5 A^ 
V MV = M V„ - d„V M + [V„, V„] + A] (85) 
A — D V A - D v A 

where we define D%f = dpf + [V M , /] 

A similar decomposition for the scalars fields is introduced: 



7T = P t (/> t (86) 



AS^ is then cast in the form: 



AS$ = 4 / Tr 



1 



d^VpA,, + V»dvVpAa + ~VpV v V p A a + -V^AvApAa 



\{D V ^ f + iA^, D^a} + Al(a 2 + tt 2 ) + 1 -{A„ tt} 2 



(87) 



3.3 Bardeen anomaly 

The theory we considered up to now is very general, so for instance one should be 
able to compute the Bardeen anomaly [Bl| . This requires a suitable modification 
of the non naive counterterm part AS^ A \ Let's consider the case in which 
G = Gr x Gl- In our model this is made by choosing that half of the generators 
and half of the generators t a L vanish in such a way that t^A = t^t b R — 0. 
Let's consider moreover the case in which the left and the right representations 
of G are the same so that the vectorial representation of G turns out to be 
defined. 

Decompose the gauge parameter as 

e = it a e a = a + 75/? (88) 

Now Vp and Ap as well as a and /?, which are respectively the vectorial and 
axial gauge parameters, can be considered as independent quantities. The gauge 
transformation 

SAp = -dpe + [e, Ap] (89) 
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decomposes into vectorial and axial gauge transformations. 
Under vectorial gauge transformations one has 



6 a V fl = -D^a ; S a A^ = [a,A p ] (90) 

and the scalar fields introduced in eq.(|S6|) transform homogeneously. 
Under axial gauge transformations one has 

S p V li = \ft,A lt ] ; 5f,Ap = -DV0 (91) 

Under a vectorial gauge transformation one gets 

S a S bare = A(a) (92) 

where A{a) is the vectorial part of the anomaly j7^ ) which indeed can be de- 
composed as 



A = A{a) + A(P) (93) 



where 



(94) 



AO =~\J ^ pa Tr[d^(d p A Rp A Ra + X -A Rv A Rp A Ro 
d v A Lp A La - \^A Lv A Lp A La 
On the other hand under an axial gauge transformation one gets 

—part) ^\P) ^Bardeen 
A B ardeen =4 J Tr {/3 [~ V ^ V pa + ^A pv A pa - ( 95 ) 

2 8 
~^{A^Ai/\^ P a ~1~ A^,\^i/ P A(j ~\~ V[ivA p A<j^ ~\~ —A p AisA p Aa 

One can then modify the renormalization of the product of the currents by 
subtracting out all the terms depending on the Levi-Civita tensor in the finite 
part of the bare action ( |S7| ) and correspondingly by introducing their gauge 
variation in eq. (|75|) . Being T unchanged the effect of subtracting eq.(|9^) is 
to recover the vectorial gauge invariance, the subtraction of eq.(^5|) put the 
anomaly in the Bardeen form |3l| . 

Observe that the e-dependent part of ( |87|) is not specific of the BMHV reg- 
ularization, being fixed its gauge variations (|92|, |95| ) . 

The e-independent part of ( ]37| ) is specific of the BMHV regularization; for 
instance it takes a different expression in [32| , where a different renormalization 
scheme is used. Using @ with V M = {A a p (t R + t L ) and A p = ^A a p (t R ~ t L ) 
and making the non-minimal choice z ab — ^S2(F) ab we find agreement with the 
finite counterterm computed in fllf using the Bonneau identities . 

3.4 Two loop 

At two loops similar computations could be performed. We have computed the 
marginal Wilsonian action in the two gauge field sector, that is the coefficients 
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(a) (b) (c) (d) 

Figure 7: Graphs of the two-point gauge field Green function. 
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Table 6: Contributions to the two-point gauge Green function: the numbers 
include the contributions of the two-loop graphs with the corresponding one- 
loop subtraction and sum over external legs permutation 



af and af of eqM 
contributions to a" 
tensors: 



at two loop. The graphs involved are shown in fig. [?]; their 
are expanded on the basis of invariant symmetric 



and a?; b 



K ? = \ ("• ['li'.te' 

K? - i ( 



,a ,b vyt 



■tr 
tr 



(« fb y\ v 



(96) 



and are shown in table ^. In order to check the first Ward identity relation 
in ( |Sl| ) we have also computed the graphs shown in fig. |s| and the results are 
collected in table ^[ Actually we did not consider in figs. 0, || the graph with 
only scalar internal lines and quartic scalar vertex: it is finite and rigid invariant, 
therefore it does not contribute to the bare action and it has to match separately 
its own T insertion. 

From the pole part of graphs in fig. [?] we computed the gauge invariant part 
of at two loops: 



S (S = j \^u[--Kt]Fl v (97) 

finding agreement with p^ ]. 

From the non-invariant part we derived the following contribution to the 
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(e) 

Figure 8: Tau-gauge graphs 
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Table 7: Contributions to the one-point T Green function. 
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bare action: 



A *(2) - 



2 53 N 

3^ - 36 ) {A^PMDlfiA^P,} + 2A,{D v v fA,P l P l - 
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-r / Tr 



|-^(^)W-(i + |](v,/Vr 



£ - £ ) 0V) 2 (Pi) 2 +* ( 7" ~ "7T ) *V&-VP 



28 29 
9e ~ 27 
16 35 
9e + 27 



8 13 
3e 9 



7 ^ Vui/^T-uy iS".; -f^i 



4 14 
9^ + 27 



1 



3e 9 
V P-S-A 



V A PS 



We used the notation of eqs. (|84| , |85| ); the r.h.s. of eq.(|98|) is not simply 
quadratic in the gauge fields because we completed it in order to have an ex- 
pression gauge-invariant under vectorial transformation. This fact has been 
possible because the derivatives d^Vv of the gauge fields V p appear only in 
the combination d[^V v ] as expected, since dimensional regularization and ours 
renormalization condition respect vectorial symmetry. 

We checked the non renormalization theorem of the anomaly only in the 
two-gluon sector of eq. (75). We have computed the part of T from which the 
anomaly might arise: 



f = / d ll e a (x)d v A b p Ale^P"T a 



:ibc 



(99) 



where T a b c is an invariant tensor symmetric in the last two indices. The non 
renormalization theorem requires that the completely symmetric part T( abc y 
which cannot be eliminated by a local counterterm, must vanish. 

Fig. ^| shows the graphs which give contributions to Ti^y Each contribution 
is decomposed on a suitable basis of symmetric tensor, we have chosen a basis 
in which 9fj never appears explicitly (symmetrization in the indices a, 6, c is 
understood): 



T { 2 = Tr[Y^Y k t R t b R t R ] Tr[Y k Y£tlt\tl\ 
Tal = Tr[Yle L Y k t R t R ] - Tr[Y k t R Ylt b L tl\ 



(100) 



The results of our calculations are summarized in table g: as expected the 
sum of the contributions of every column vanishes, so that the Adler-Bardeen 
theorem is verified. 

One could address the question whether T abc and not only its symmetric 
part is actually vanishing. For G = SU(N), using the Young tableaux one can 
easily prove that invariant tensors X abc with mixed symmetry do not exist; ^ 
for more general groups see f33| l . 

3 We thank C. Destri for explaining this point to us. 
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(b) 



(c) 



(d) 



(e) 




(f) 



(g) 



(h) 



Figure 9: Two- loop graphs contributing to the anomaly 



T. 



(i) 



abc 



T 



abc 



graphs 



6797 ■ _ 4688 ■ 
21870 6561 



18119 ^ _ 12880 iv 



21870 



6561 



254 ■ , 2080 ■ 
'243 ^ 729 lu 



135 243 v 



1861 
1215 



1192 , 
729 



. 55441 j - J§-i v 

204120 2187 



10931 „■ _ 308 ■ 
204120 2187 lv 



135 ~ 243 



1861 ■ _ 1192 ■ 
1215 729 lu 



Fig. g 
Fig. gb| 
Fig. gSJ 
Fig. gd)| 



Fig. ^ 

Fig. |m| 



Fig. |9(g_ 



Fig. £(hj 



Table 8: Two- loop coefficients for the anomaly 
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4 Concluding remarks 



We have shown that Wilsonian methods are useful in the BMHV dimensional 
renormalization of theories with chiral symmetries; the two-loop renormalization 
of the most general Yukawa theory becomes straightforward in this scheme, 
while in a more conventional approach, based on the verification of usual Ward 
identities, it is a non trivial task. 

In our approach one renormalizes the effective Wilsonian action along a 
flow on which the subtraction integrals are easily computed using standard 
techniques. 

There is some arbitrariness in the choice of this flow; we worked with the 
n = 2 flow, which is the simplest one from a computational viewpoint: it is 
quite obvious that calculations with the n > 3 flows are more complicated, 
although the bare action is the same. On the other hand the n = 1 flow, which 
coincides with the auxiliary mass method, does not respect renormalizability so 
the procedure described in this paper fails down. It is a convenient method in 
nonchiral theories where minimal subtraction respects the Ward identities. 

In the future we intend to apply our formalism to chiral gauge theories. 
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Appendix. Two-loop bare Yukawa action 



We have written the naive part of the counterterms in the Tables; here we report 
the remaining counterterms, necessary to recover the rigid chiral symmetry in 
the BMHV scheme. 

The two-loop scalar quadratic counterterm is 



Ac«(p) = ctp 2 + c™ 2 ' kl (n)m 2 kl +f terms 



c?. = Tr 

n 



2S l P k S J P k +(-- + - PiS k PjS k + ( — - — ) (SiSkPjPk + 



StPkPjSk + PiPkPjPk) + ( 1 - - ) SiSjPkPk + 



c 

7 4 \ f 13 16 

g - Y e ) ( s > s kPkPj + s l p p k s k + PiPjSjSk) + ( — - — 



= Tr[(3 _ i ){SjSkPiPl + SjPkPiSt + PjSkSiPt + PjPkStSt) - 
J e 

SPjSkPiS! + - ^PjPkPiPt] + - ^PjPkPtPi + 

(2 - *)(SjS k PtPi + SjPkPiSi + PjSkStPi + PjPkSiSi) + - l)P m ^]hi jmk 



and the symmetrization in i, j is understood. 

The two-loop fermionic quadratic counterterm is 

Ac(p) = ^ dj, + y> — terms] 



The two-loop fermion-fermion-scalar counterterm is 
7 



Aci = ^2 c i( n ) 



71=1 

c i (l)=iy k [(- 7 -+*)Tr(P j P k )S i + 

12 22 32 

^5(2 - -)Tr(5 j 5 fe )P i +^75(2^ - -)Tr(P,P fe )P] yj 

ci(2) = iy k [h 5 (\ - \)(SjSiPj + PjS^) + 

i lb 2 -S 3 P t S 3 + {-\ + ^(SjPiPj + PjPiSj) + i lh { 2 - - ^PjPiPAvk 

Ci(3) - iy*^ - ^(SiYiPj + PiYjSj + 
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Ci(4) = iy k [i^{\ - \)(S t Y 3 P k + PjYkSt) - z 75 (i + ^{PiYjSk + S.jY k P t ) 

-(PiSjPk + P S k P t ) - i- /5 (l + -)(PiPjP k + PjPkPi)]yj 

Ci(5) = iykliTspSjPiSk - \s 3Vl P k - ^ > P J P l P k - ^PjyiS k )] yj 
c<(6) = 

c *( 7 ) = --^VjPkyihijki 



The two-loop quartic scalar counterterm is 

9 



n=l 

c(l) - 16 TrF m [z(^ - ^)5F m P 3 + - jj)Sy m (SyP + PFS) + 
(- + H)py m p3 + py-^f (gytp + p3 + py S y + 

- i)pr m (5rp + pys) + »(^ - Jj)py ro (sys - pyp) + P . c . 

c(2) = 4 Try ro [^(ytyy r tyrt _ y2y m y2) _ 

8P 2 y m (y 2 + ^p 2 ) + (syty - p 2 )yIp 2 \ + p.c. 

8 

C (3) - rry m [(- - 3)(y 4 y m - (yty)^) + i2P 4 y m - (- + ^)p 4 y,t + 

6 6 O 

oo o 

(— - -){2iSYP 2 P m + 2iPS 2 PP m - P 2 Y^SYl + P 2 YSY m )\ + p.c. 

6 o 

AO 

c(4) - Tr[-(— - 12)(S ro SPP„ + S m P 2 S„ + 

P m S 2 P n + P m PSS n ) — ( — ) P m P 2 P n ]h mn 

c(5) = Tr[-(— - ±|V m PP n P - 48S m PS„P - 
24 

(— — 18)(S m SP n P + S m PP n S + P m S S n P + P m PS n S)]h mn 
6 7 

c(6) = Tr[(- - -)P m P„]/i mp /i„ p 
c(7) = c(8) = c(9) = 



where the indices i,j, k, I, which are totally symmetrized, are understood. 
p.c. represents the terms obtained with the substitution rule: P — > — P, Y — > yt , 

yt -» y. 
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